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q \ Abstract 

CN ' The phenomenological investigation of noncommutative space-time in the laboratory 

frame are presented. We formulate the apparent time variation of noncommutativity 
parameter 9^ in the laboratory frame due to the earth's rotation. Furthermore, 
in the noncommutative QED, we discuss how to probe the electric-like component 
Ge = (#oi) $02, #03) by the process e~e + — > 77 at future e~e + linear collider. We may 
determine the magnitude and the direction of 9e by detailed study of the apparent 
time variation of total cross section. In case of us observing no signal, the upper 
limit on the magnitude of 9 E can be determined independently of its direction. 

(N 
O . 

{Sj , 1 Introduction 

o 

The early study of noncommutative space-time was presented by Snyder]^] in 1947, with respect to 
the need to regularize the divergence of quantum field theory. In Snyder's work, it was suggested 
that the divergence may be regularized by an elementary unit of length induced by the noncommu- 
tativity of space-time. Snyder's basic idea was the extension of the quantization of phase space in 
quantum mechanics. Furthermore, noncommutativity of space-time may arise from string theory 
in the specific low energy limit Q. The noncommutative space-time is characterized by operators 
5_j ■ X^ satisfying the commutation relation 



(N 



[X fl ,X v ]=i9 llv (1.1) 

where 9^ v is antisymmetric constant matrix, 9^ v = —9 vtl and \X p ,9uJ\ = 0. And 9^ u have dimen- 
sion of (Length) 2 . Therefore, ( |1.1| ) introduce the elementary unit of length in the theory, such 
as Planck constant in quantum mechanics. Nonzero constant matrix 9^ v may violates Lorentz 
invariance. Lorentz violation in noncommutative quantum field theory have been studied |J. 

It is known that QED in noncommutative space time (NCQED) [Q] is invariant under the 
noncommutative version of U(l) gauge transformation and is renormalizable at one loop level 0, 
|, §. Axial anomaly §] and CPT invariance J?]] in NCQED have also been studied. There are 
several phenomenological study on NCQED for low energy experiments (§], |j| [l(| 11]. Assuming 



9^ is constant in the laboratory frame, a lower bound on noncommutativity scale A^vc have been 
found to be A^vc > 100GeV|lO|1 in order that the result of Lamb shift is consistent with the 
ordinary quantum mechanics. Other limit on noncommutativity parameter have been found to 
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be 9 < (10TeV)~ 2 , if 9^ v = 9e^ v , by an analysis of noncommutative Aharonov-Bohm effect |fLl|| . 
High energy phenomenology in NCQED has also been studied for several processes at future linear 
colliders [p"2|, |l3|| . Moreover, phenomenology relevant to Standard Model (SM) like interactions in 
noncommutative space-time have also been studied [14, 15] on the assumption that we may obtain 



SM-like interaction in noncommutative space-time by usual procedure replacing every products 
of fields with the star product. In those previous studies, however, the direction of 9^ v have 
been assumed to be fixed to the laboratory frame. Such an assumption might be justified, if 
measurements would be given by the data set suitably averaged over time and also over polar 
angle distributions. 

The 9fj, u , however, may be considered as an elementary constant in the nature. And there may 
exist a class of specific coordinate system in which the direction of 9^ u is fixed. It is likely that 
such a coordinate system is fixed to the celestial sphere. 

On the contrary, the laboratory frame is located on the earth and is moving by the earth's 



rotation. Therefore, as was mentioned injL^, 15, 16], we should take into account the apparent time 
variation of 9^ v in the laboratory frame due to the earth's rotation when we discuss phenomenology 
for any experiment on the earth. In this paper, we will consider the effect of apparent time variation 
of 9^ v in the experiments due to the earth's rotation seriously. 

If an anisotropy due to noncommutativity of space-time exists, probing the specific direction 
of 9^ v and measuring the magnitude of elementary unit of length are very interesting tasks from 
both experimental and theoretical aspects. We may determine the direction of 9^ by the analysis 
taking into account effects of time variation of the measurements. 

This paper is organized as follows. In section 2, we present the parameterization for 9^ v 
including the effect of the earth's rotation. In section 3, we make some comments on the time 
dependent cross section and we define the time averaged cross section. In section 4, we show several 
numerical results and we discuss how to prob 9^ v at future linear collider experiments. Finally, we 
conclude our result and discussion. 

2 Expression of 9^ v in the laboratory frame 

The noncommutativity parameter 9^ can be classified into two parts. One is the electric-like 
component 9e = (#oi> #02, #03)- Another is the magnetic- like component 9b = (#23; %L> $12)- Those 
elements can be determined when a coordinate system is chosen. In the specific coordinate system, 
both 9e and 9b should be constant vectors. Hereafter we call such a coordinate system a "primary" 
coordinate system. It is feasible that we take a set of coordinates fixed to the rest frame of the 
cosmic microwave background(CMB) as a "primary" coordinate system. According to COBE 



experiment [17], the boost of the solar system for the CMB rest frame is about 370km/s. This is 
about 0.12% of the speed of light in vaccum. Moreover the speed of the earth in solar system is 
about 29.78km/s. Therefore the effect of the boost to the measurement of 9e and 9b are small 
enough to neglect in comparison with the detector resolution in the collider experiments. And 
we may consider that the CMB rest frame is fixed to the celestial sphere approximately. Thus, 
hereafter, we assume that a primary coordinate system and also each direction of 9e and 9b are 
fixed to the celestial sphere effectively. 

At first, we introduce a primary coordinate system which is Cartesian coordinate system. The 
Z axis is along the axis of the earth's rotation and the positive direction of Z axis points to the 
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Figure 1: "primary" coordinate system (X- 
Y-Z). The axis X point to the vernal equinox 
Tj2ooo.o- The electric- like component Oe of 
Ofxv is also shown. The direction of Oe is pa- 
rameterized by constant angle parameters rj 
and £. 
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Figure 2: Arrangement of laboratory coor- 
dinate system (x-y-z) for an experiment on 
the earth in the "primary" coordinate system 
(X-Y-Z). 5, a and u are constants. 



north. The axis pointed to the vernal equinox (Tj2ooo.o) is labeled X. We take X-Y-Z system as 
the right-handed system. Figure [l] shows the sketches of the primary coordinate system and the 
direction 9e parametrized by 77 and £. 

Let ex, &y and e~z be the orthonormal basis of the primary coordinate system (X-Y-Z). Then 

Oe = Oe (ex sinr/cos^ + e~y sinr/sin^ + ezcosrj) , (2-1) 

where < 77 < 7r, < £ < 2-7T and 6 E = \0 E \- To be exact, this X-Y-Z coordinate system moves 
slightly owing to the earth's precession. Since the period of the earth's precession is about 2.6 x 10 4 
years, the vernal equinox is moving by about 0.014 degree/year. Therefore, we can neglect the 
earth's precession during the term of experiments. 

On the other hand, the usual coordinate system for experiments is fixed to the detector. We 
label each axis of such a coordinate system by small letter (x,y,z). As an example we consider 
an e~e + collider experiment. The origin is set at the interaction point. The z axis is along the 
direction of e~ beam. The horizontal and vertical axes are labeled x and y respectively. The 
x-y-z system should be the right-handed system. Hereafter we call this coordinate system the 
"laboratory" coordinate system. 

As is shown in figure [2], we parametrize the location of an e~e + experiment on the earth by a 
latitude 5 of the detector site, the angle a between direction of z axis and the meridian at detector 
site, and the angle £ between X-Z plane and y-Z plane. The angle a is measured counterclockwise 
from the north. [ 

1 Our definition of the angle a is opposite to the definition of usual azimuth in astronomy. We define angle a as 
it increases with a positive rotation in the right-handed system. We may call the angle a the counter-azimuth. 
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Let and be the orthonormal basis of the laboratory coordinate system (x-y-z). The 

transformation between (ex, ey, e~z) and (ej, e^, ej) is given by 





(2.2) 
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(2.3) 



with — 7r/2 < 5 < 7r/2 and < a < 2ir, where we use the usual abbreviation, = cos£, etc. 
Hereafter we take the orthonormal basis of the laboratory coordinate system as the usual way, 
= (1,0, 0) T , = (0, 1,0) T and el = (0,0, 1) T . Then, in the laboratory coordinate system, the 
orthonormal basis of the primary coordinate system can be written as 

/ -CaCQ + SsSaSQ 
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(2.4) 



Note that in the laboratory coordinate system the direction of Z axis, namely the axis of the 
earth's rotation, is given only by the location of e~e + experiment (<5, a). For example, (<5, a) 
of LEP experiments @ are approximately (46.15°, 40°) for OPAL, (46.15°, 130°) for ALEPH, 
(46.15°, 220°) for L3 and (46.15°, 310°) for DELPHI. Therefore the arrangement of e% in the 
laboratory coordinate system of the LEP experiments are different each other. 

The angle £ increases with time t owing to the earth's rotation. A detector site will return to 



the same direction by a sidereal day, T^ ay =23h56m4. 09053s [19]. Therefore, we may take 

C = UJt With LU = 2lT /Tday, 



(2.5) 



by setting t = when the detector site is on the X-Z half plane with X > 0. 

Substituting (|2.4| ) into fl2.1|) , we find the expression of 9 E in the laboratory coordinate system, 
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(2.6) 
(2.7) 

(2.8) 



where 9es is the projection of 9e onto the Z axis and is the stable part of E in the laboratory 
coordinate system. Oev is the time variation part of 9e- The direction of 9ev revolves on the 9es 
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Figure 3: Two typical time variation of Oe in the laboratory frame, (a) for 77 < ®Lab and 
(b) for rj > <d Lab , where cos0 La6 = c s c a . 



by a period T^ay This is the apparent time variation due to the earth's rotation. Angle parameter 
£ appears in the expression of Oev as the initial phase for time evolution. It is easy to show that 



Oes 



Oes ■ &EV 



cos 77, 
= 0, 



Oev 
► 2 

Oes + Oev 



sinij, 



0e 



(2.9) 



Therefore the magnitude of each vector QeSi Oev and 6e is independent of time. 

Let ® Lab be the polar angle of Oes in the laboratory coordinate system. From fl2.8| ), we find 
cosQ Lab = c s c a . We may classify the apparent time variation into two typical cases, 77 < @ Lab and 
77 > Qi ab - Figure ||] shows the two cases for the apparent time variation of Qe in the laboratory 
coordinate system. Let <3? E be the azimuthal angle of Qe in the laboratory coordinate system. In 
the case (a) 77 < Q Lab , &e varies within the region of (<3?£ ax — $™ n ) < ir. On the other hand, in 
the case (b) 77 > ®Eabi varies within the whole region. Therefore we may expect that some 
typical differences exist between the case (a) and (b) in the angular distribution for a process, for 
example e~e + — * 77, which are affected by the space-time noncommutativity. 

The magnetic-like component Ob is also parametrized by the same way. In general, however, 
both the direction and the magnitude of Ob are different from those of 0e- Therefore 0^ can be 
parameterized by six parameters, four angles and two magnitudes of Oe and Ob, e = \6e\ and 
9 B = \9b\- In the primary coordinate system, 



1e = Oe (ex sin rj E cos £ E + ey sin rj E sin £ E + e~z cos rj E ) , 
9b = Ob (ex sin rj B cos £b + ey sin tjb sin £ B + ez cos t)b) ■ 



(2.10) 
(2.11) 



By using (2.4), we can obtain the expression of Oe and Ob in the laboratory coordinate system. 
We should take Oe and Ob as model parameters rather than the energy scale Ae = l/y/0~E and 
Ab = 1/V0~b~. 



e e ' 



77 in NCQED 



A field theory in noncommutative space-time can be described equivalently by a field theory with 
commutative space-time variables in which every products of fields are replaced with the star 
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product of fields. The star product is defined by 

(3.1) 



y=z=x 



f*g(x) = exp r-d^e^ f(y)g(z) 

where x, y and z are ordinary commutative variables. 
NCQED action|]] is given by 

S = J d A x (- ^F^ * + * — m^f^\ (3.2) 

where F^ u = d^A v - d u - ie (A^ * A u - A u * A^). And covariant derivative for the matter fields 
is given by = — ieA^-kfy. We need nonlinear terms in field strength F^ v to keep NCQED 
action invariant under noncommutative U(l)± gauge transformation, 



Ap - A' II = U(x)*A^U~\x)- l -(d ll U(x))*U' 1 (x), (3.3) 

tf(x) -> V'(x) = U(x)*W(x), (3.4) 
V(x) -> ¥(x) = ^(x)-kU~ 1 (x) (3.5) 



where 



U(x) = exp(ia(s))* = £ M^M! (3.6) 

n=0 n - 

and f/(a;) * f7 _1 (») = U~ 1 (x) * C/(x) = 1. The NCQED action (|J) is invariant under the 17(1)* 
gauge transformation. 

In NCQED, we consider the pair annihilation process e~ (k 1 )e + (k 2 ) — > 7(^1)7(^2) a t future 
e~e + linear colliders. Each momentum is taken to be 

^ = ^(1,0,0,1), ^ = ^(1,0,0,-1), # = j£ = (^,-?). ( 3 - 7 ) 

where = (-^/s/2)(s e c ( j ) , s^s^, c e ). and <^> are polar angle and azimuthal angle of final state photon 
in the laboratory coordinate system. The differential cross section for e~ (fc 1 )e + (k 2 ) — > 7(^1)7(^2) 
in the center of mass system is given by 



da a 2 ( t u t 2 + v? 

= — — 1 4 

dcos Odd) 4s \ u t s 2 



Ui V Lb O . 1 /> 

+ - - 4 ^ sin 2 A wc I , (3.8) 

Anc = rffeS = _^^ (39) 

where s,t and it are usual Mandelstam variables, s = (k 1 + k 2 ) 2 , t = (k 1 — Pi) 2 and u = (k-> — p 2 ) 2 - 
The E is given in (|2.6| ), fl2.7|) and fl2.8|) . When A^vc = 0j differential cross section ( |3.8D corresponds 
to the differential cross section in QED. 

Since two photons in the final state are identical, we cannot distinguish two configuration (0, dS) 
and (tt — 0,tt + dj). Therefore, we must get the sum of the differential cross section for (0,d)) and 
(ir — 0, ir + (p). Then the observable is 

d(J ° bs ~ ^ "^ + _^(7r-«U + tf. (3.10) 



d cos Odd) d cos Odd) dcos 
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Note that (|3.10| ) is defined in the region < cos 9 < 1 and < cf> < 2-ir. It is easy to show from 
and ( |3.9| ) that A^ci^ — 9 ,ir + <j)) = —A^cifi,4>) and 

' h -(^nc) = -n^i-A-^Nc)- (3-11) 



d cos 9d(fi dcosi 

Moreover, this imply that the differential cross section ( |3,8|) and ( |3,10| ) are symmetric for the change 
of the sign of 9 E , 9 E <-> —9 E . Therefore we cannot distinguish between (r/, £) and (tt — 77, tt + £) by 
observing the process e~e + — ► 77. There is two-fold ambiguity for the determination of (r/, £). 

We can see from (|3.8| ) that NCQED effect to the differential cross section of e~e + — > 77 always 
gives the negative contribution, moreover, from fl3.9|) we find 



(AatcI = max(Ajvc') = — 0# if II &E, ^ ^ 

A NC = if -f JL 6^. 

This means that, when we compare NCQED prediction with QED prediction, the deficit of the 
differential cross section appears around the specific direction in which J? is almost parallel to 
9 E - Furthermore, such a specific direction varies with time in the laboratory coordinate system, 
as we have discussed in previous section. Therefore, in general, observables for e~e + — > 77 in the 
laboratory coordinate system have time dependence even for the total cross section. 

We may consider that the measured value for observable by collider experiments is usually 
given as a mean value. And such a mean value should be compared with NCQED prediction 
averaged over time. Taking into consideration that the period of time variation of the observables 
in NCQED is the sidereal day Td ay , we introduce the time averaged observables as follows; 

da \ 1 [ T day da„ 



d COS 9d(j) I T Tday Jo d cos 



obs 



dt, (3.13) 



ih X ' / '' : '" da °>° dt, (3.14) 



d cos 9 1 T T^y Jo d cos 9 
da\ 



> / T Tday JO 

\ fTday 



1 ' "^11. (3.15) 



{<t)t = 7f~ / Vobsdt, (3.16) 

-i day JO 



where 



d(T obs _ f 2 *^ d(T obs 



(3.17) 



d cos 9 Jo d cos 

d(p Jo d cos 9dq> 

aa obs 



a obs ^ I d(cos9) I dtj^fjj. (3.19) 





The polar angle cut is denoted by e (0 < e < 1). It is easy to see that (<t)t = a obs which is usual 
total cross section, when a obs is independent of time. 

Note that we have integrated out the £ dependence of the observables by taking average over 
time, since £ play a role of initial phase for time evolution. Therefore 9 E and angle 77 may be 
determined by the time averaged observables. 
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Figure 4: Time averaged azimuthal angle distribution for rj = 0, 7r/6, it/ 3, it/ 2, 2tt/3, 
5tt/6, it. We set the laboratory coordinate system by taking (5, a) = (vr/4, vr/4). We take 
y/s = 500GeV and 6 E = (500GeV)~ 2 . 



4 Numerical Results 

We show several characteristic results in NCQED and also discuss how to prob by using ob- 
servables in the laboratory coordinate system. We set the laboratory coordinate system by taking 
(8, a)=(-7r/4, 7r/4). The cut for cos# is taken e = 0.2. 



4.1 Azimuthal angle distribution 

Anisotropy of azimuthal angle distribution of e~e + — > 77 is predicted in NCQED even if we 
consider the time averaged distribution (da/d(j)) T . Figure || shows (da/d(p) T for 9 E = (500GeV)~ 2 
and several values of rj. We take -y/s = 500GeV. 

We see from figure f|that the curves of (da/d(p} T are sensitive to the value of rj around <j) — 0-7-7T 
and also almost independent of the value of r] around <p ~ 1.7-7T. Furthermore {da/d(f)) T is almost 
flat around <fi ~ 1.7tt for any rj. 

Those specific angles 0.77T and 1.7-7T can be interpreted as the azimuthal angle of 9 ES and —0 ES 
in the laboratory coordinate system. In other words, those are the azimuthal angles of north pole 
and south pole of the celestial sphere. Since we take (5, a) = (7r/4, vr/4), the azimuthal angle 4> ES 



of 9 ES can be derived from (pll) as follows 



^-^S-T. • ^-JT=S-^ (41) 



-5 u a v V 5 a 

then (f> ES ~ 0.77T and the azimuthal angle of —0 ES is given by (f> ES + tt ~ 1.77T. 

We also see from figure || that each input rj and it — rj gives the same distribution of {da/d(j)} T . 
This is because the differential cross section of e~e + — ► 77 is symmetric for 9 E <-> —0 E . 
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Figure 5: Apparent time variation of total cross section for rj = 0, vr/6, tt/3,tt/2, 2tt/3, 5tt/6 
and 7T. We take (5, a) = (vr/4, tt/4), y/s = 500GeV and # E = (500GeV)- 2 . 



We may determine 77, except for two- fold ambiguity for rj and ir — r], by fitting the shape of curve 
of (da/d(f)) T , especially around <j) ~ 4>^ s = 0.7tt. Also we may determine 6 E almost independently 
of 7} by measuring the deficit of {da/d(j)) T compared with QED prediction around <j) ~ 1.7n. 

4.2 Time dependent total cross section 

In order to determine £, we need to trace the apparent time variation of observables due to the 
earth's rotation. Since total cross section a obs depends on 6 E , rj and also £, we may expect that 9 E , 
r\ and £ could be determined by measuring time variation of a obs precisely, except for the two-fold 
ambiguity for (rj, £). 

Figure [| shows a obs as a function of cot — £ for 9 E = (500GeV) -2 and for several values of rj. 
We can see from figure || that the shape of curve is sensitive to the r/. If the time variation of total 
cross section is observed, we could determine both magnitude and direction of 6 E by fitting the 
NCQED prediction of a obs with the data in three parameters space (0 E ,r],£). The magnitude 9 E 
and the angle r\ may be determined by the fitting of both the magnitude and the shape of curves 
of a obs . The £ may be determined by the measurement of the phase of time evolution of cr obs . 
However, as is mentioned previously, we cannot distinguish (77, £) from (7r — 77, 7r + £). For example, 
the graph of a obs for 7/ = ir/3 is identical to that for r/ = 2ir/3 shifted the phase £ to £ + ir. 

Although we may determine 9e by tracing the time variation of the differential cross section 
of e~e + — > 77 instead of total cross section, we can easily imagine that such an experiment needs 
very large luminosity. This is because we must divide not only the phase space but also the time 
distribution into many bins, in order to trace the time variation. Therefore, in the determination of 
9 E , ij and £, we had better probe the time variation of total cross section in the early experiments 
at e~e + linear colliders. 
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QED 

6 E =(1000GeV)" 

6 E =(600GeV)" 2 



9 E =(500 GeV) 

6 E =(400GeV)" 2 
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Figure 6: Time averaged total cross section of e e + — ► 77 for ^ = (400GeV) 2 , 
(500GeV)~ 2 , (600GeV) -2 and (lOOOGeV)" 2 are shown. Horizontal axis is taken to be 
angle r\. QED result is also shown. We take (5, a) = (ir/A, 7r/4) and = 500GeV. 



4.3 (da / d<fi) T vs. a o6s 

We can see from figures ^| and || that (da/d<f>) T and cr f, s show opposite behavior for each input 
value of 77. For example, when 77 = or tt, we may observe large variation of azimuthal angle 
distribution (da/d(j)) T . In this case, we find no time variation of <J bs- On the contrary, when 
77 = 7r/3 or 27r/3, since the variation of (do /d(j)) T is very small, we may observe the flat distribution 
in the experiments. In this case, we find large time variation of cr 6s- Therefore we may expect 
that non-uniform distribution due to NCQED effect should appear in the (da/d<p) T and/or a b s 
for any value of 77, if 9 E is large enough. 

4.4 Time averaged total cross section 

Finally, we consider what we can measure by the time averaged total cross section (ct)t- Figure^ 
shows ((j)t as a function of 77. It is easy to see that (c) T is almost independent of 7/. In case of us 
observing some deficit of (cr) T , we may determine 9 E independently of 77 by measuring (a) T . 

On the other hand, in case of us observing no signal, we may obtain the upper limit on 6 E . 
The lcr deviation for total cross section in QED, ctqed, can be estimated by J '(Jqed/ ' L. For 
\fs = 500GeV and e = 0.2, we have (28.89/\/i) where L is the luminosity given in fb — 1 . In this 
case, an expected 95%CL upper limit on 9 E is found to be 

9 E < (600GeV)- 2 for L=100fb- 1 . (4.2) 

Furthermore, since \oqed — (o)t\ ^ ( s ®e) 2 for \s9e\ < 1, we may estimate 95%CL upper limit on 
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9 E for arbitrary L from (4.2) as follows; 




lOOfp- 
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e E < (eoocev) 



(4.3) 



For example, we find E < (800GeV)~ 2 when L=1000fb" 1 . 
5 Conclusion and Remarks 

We have presented phenomenological formulation of the apparent time variation of noncommu- 
tativity parameter 9^ v in the laboratory coordinate system. In our framework, the laboratory 
coordinate system have been taken to be a familiar coordinates to the collider experiments, and 
the primary coordinate system fixed to the celestial sphere have been introduced. We have shown 
the transformation formula between the primary and the laboratory coordinate system, and also 
shown the expression of 9 E in the laboratory coordinate system. The formulation presented in 
this paper is applicable to the study of models which predict an intrinsic direction of the space- 



As an example, we have applied our formulation to NCQED and discussed the determination 



parameterized by 9 E , rj and £ in the primary coordinate system. We have shown that 9 E may 
be determined by the detailed study of the time dependent total cross section, though two-fold 
ambiguity in the parameter space (77, £) remains. To determine £, we need to probe the phase of 
the time evolution of cr & s . I n case of us observing no signal, probably this is the most realistic 
case, we may obtain the upper limit on 9 E independently of the direction of 9 E . 

So far we have considered one experiment with (<5, o)=(7r/4, vr/4). If there are several detector 
sites in the e~e + collider experiment and the direction of e~ beam in each site is set to be along to 
the different direction, such as four LEP experiments, then the angular distributions of e~e + — > 77 
and the time variation of observables should behave differently in each experiment. This is because 
the direction of 9 ES in the laboratory coordinate system at one detector site differ from that at other 
detector sites. Therefore we can expect that the combined analysis of several experiments with the 
different (5, a) play an important role in the attempt to probe the space-time noncommutativity. 

Finally, we would like to make some comments on the determination of the magnetic-like 
component 9 B . Since the process e~e + — > 77 is independent of 9 B , to determine 9 B , we must 
consider other processes relevant to 9 B , for example e~7 — ► e~7 process which depend on both 9 E 
and 9 B . The process 77 — > 77 may also be available to determine 9 B . By combining the results 
from those processes, we may determine 9 E and 9 B . We postpone the study of this matter to the 
future studies. 
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